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We consider the recent calculation Q] of the graviton propagator in the spinfoam formalism. Within 
the 3d toy model introduced in Q, we test how the spinfoam formalism can be used to construct 
the perturbative expansion of graviton amplitudes. Although the 3d graviton is a pure gauge, one 
can choose to work in a gauge where it is not zero and thus reproduce the structure of the 4d 
perturbative calculations. We compute explicitly the next to leading and next to next to leading 
orders, corresponding to one-loop and two-loop corrections. We show that while the first arises 
entirely from the expansion of the Regge action around the flat background, the latter receives 
contributions from the microscopic, non Regge-like, quantum geometry. Surprisingly, this new 
contribution reduces the magnitude of the next to next to leading order. It thus appears that the 
spinfoam formalism is likely to substantially modify the conventional perturbative expansion at 
higher orders. 

This result supports the interest in this approach. We then address a number of open issues in 
the rest of the paper. First, we discuss the boundary state ansatz, which is a key ingredient in the 
whole construction. We propose a way to enhance the ansatz in order to make the edge lengths and 
dihedral angles conjugate variables in a mathematically well-defined way. Second, we show that the 
leading order is stable against different choices of the face weights of the spinfoam model; the next 
to leading order, on the other hand, is changed in a simple way, and we show that the topological 
face weight minimizes it. Finally, we extend the leading order result to the case of a regular, but 
not equilateral, tetrahedron. 



I. INTRODUCTION 

The spinfoam formalism Q is a candidate covariant approach to a non-perturbative quantisation of General Relativ- 
ity (GR). At present, it provides a consistent background independent theory at the Planck scale, where it describes 
spacetime as a discrete quantum geometry. However, the large scale behaviour of the theory is less understood. 
Indeed, the formalism lacks a well-defined procedure to study the semiclassical limit, to define particle scattering 
amplitudes and to reproduce low-energy physics. In particular, considering the pure gravity case, we would expect to 
recover in the low-energy regime the conventional perturbative expansion described in term of gravitons. Recently, 
an important step towards this semiclassical limit has been taken with a proposal for the construction of the 2-point 
function of quantum gravity from the spinfoam amplitudes Q . It relies on the use of the propagation kernel (see for 
instance 1^1) and on a boundary state peaked around flat geometry. The proposal has been considered in 4d Q, Q 
and 3d Q Riemannian GR without matter. In both cases the large scale limit of the linearised 2-point function has 
been shown to reproduce the expected 1/p 2 behaviour of the free graviton, thus providing a first piece of evidence 
that spinfoams might correctly lead back to general relativity in the semiclassical limit. 

One of the appealing aspects of the proposal is that it can be used to compute not only the free propagator, but also 
the quantum corrections due to the intrinsic non-linearity of the theory. Such corrections have been studied in the 
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conventional (non-renormalisable) perturbative approach to quantum gravity, where they are understood as graviton 
self-energies. In the (physically interesting) 4d case, these give rise to corrections to the Newton potential (see for 
instance 0, 01 )■ However, the perturbative approach can only be considered as an effective field theory, and not as 
a complete theory, precisely because of its non-renormalisability: an infinite counterterm arises at two loops On 
the other hand, it has been argued that the cure to this non-renormalisabily could lie in the use of non-perturbative 
approaches, such as the spinfoam formalism. It is thus very interesting to compute the quantum corrections to the 
free propagator within the spinfoam formalism, both to give stronger support to the proposal of Q , and to understand 
how full quantum gravity modifies the usual perturbative expansion of gravitons. 

As a warm up exercise for the 4d case, we consider in this paper 3d quantum gravity and study the emergence of 
higher order corrections to the free propagator. As is well known, 3d general relativity (GR) has no local degrees of 
freedom. Due to this peculiarity of the 3d case, the following clarification is necessary. If the metric is quantised as a 
whole, the theory is exactly solvable and this is the basis of the topological Ponzano-Regge (PR) model that 

we consider here. However, if one treats the metric with the conventional perturbative expansion g^ v = r]^ v + h nv , 
the field /i M „ is a pure gauge quantity: consequently, the quantum theory does not have a propagating graviton [121 ]. 
It is nonetheless useful to consider the above expansion in a gauge where is not zero. Then the 2-point function, 

W^ p(T {x,y) = (0|T{V,(x)V(y)}|0), 

can be evaluated in a perturbative expansion in £p. The leading order, corresponding to the free theory, goes as 1/p 2 
in momentum space, which in 3d means a \ jl dependence on the spacetime distance; analogously one can compute 
the self-energies, and the perturbative expansion is not rcnormalisable [T2I Il3j. as in 4d. Therefore, it makes sense 
to study in 3d how the spinfoam-defined 2-point function can cure the non-renormalisability of GR by producing a 
different perturbative expansion than the usual one. 

To that end, we consider here a simple toy model introduced in 2] : the contribution to the graviton propagator by 
a single tetrahedron belonging to a discretisation of spacetime. A key advantage of this toy model is that it allows 
efficient numerical simulations, which are substantially more difficult in the 4d case. Indeed, it was the first numeric 
results of Q that gave hints of deviations at short scales from the inverse power law There is a precise reason 
to expect that such corrections would show a different structure than the conventional perturbative expansion, and 
that is the modification of microscopic structure by quantum geometry. To see why this can be expected, recall first 
that since we are dealing with 3d Riemannian GR, we consider the Ponzano-Regge spinfoam model. It is constructed 
in terms of half-integers (spins) labelling SU(2) irreducible representations (irreps), variables which are related to 
geometric quantities. The large scale limit then corresponds to the large spin limit. In that regime, the model is 
dominated by the path integral for Regge calculus, with a specific choice of measure. Regge calculus is a discrete 
approximation to classical GR, which captures the non-linearity of the theory. Therefore, in this limit, one can 
still expand around a flat background, and consider only small perturbations around it. The leading order, which 
reproduces the \ jl free graviton, is obtained by keeping only the quadratic term in the Regge action, and the trivial 
background measure. Higher orders in the action and in the measure give corrections with the structure of the usual 
continuum perturbative expansion. 

Thus far, one would conclude that the spinfoam formalism simply produces a somewhat discrete version of the 
conventional theory. However, because the Regge path integral only emerges in the large spin limit, the full theory 
predicts further sources of corrections, arising from the exact, non Regge-like, quantum geometry. Heuristically, these 
corrections are remnants at large scales of the discrete microscopic structure of the spinfoam geometry. These other 
corrections have no analogue in the conventional theory, and their presence is the reason why we expect spinfoams to 
introduce new features in the perturbative expansion. Indeed, it has often been suggested that the microscopic picture 
of quantum geometry emerging from spinfoams would affect in some way the large scale dynamics, possibly giving 
better finiteness properties. Here we show a clear example of how this could happen. In particular, it is interesting 
to study the interplay between the two different kinds of corrections. We compute the next to leading (NLO) and the 
next to next to leading (NNLO) orders, and we show that the quantum geometry corrections only enter the NNLO, 
and moreover reduces its magnitude. This is the main result of this paper. 

A key techical input to the expansion is the choice of boundary state entering the definition of the 2-point function. 
In principle, the boundary state should be given by the vacuum state of the theory, but in the absence of a well-defined 
prescription for it, we must make an ansatz. In [l| a Gaussian peaked around the (discrete quantities representing) 
the intrinsic and extrinsic geometry was suggested for the boundary state, however there has been little investigation 
of alternatives. So far, the only requirement on the boundary state is that it be a good semiclassical state; namely, 
that the relative uncertainties of the geometric quantities vanish in the large spin limit. As the structure of the state 
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is the same in both 3d and 4d, it is useful to investigate its properties in the toy model here considered. The simple 
Gaussian introduced in [ij has two problematic aspects: first, the variables representing the intrinsic and the extrinsic 
geometry are not conjugate to each other in a mathematically well defined way; second, it has an intricate (D-structure, 
due to the use of a complex phase. Here we show how both issues can be solved, by working with a smooth function 
on SU(2), and using the harmonic analysis over the group. We propose a new state and show that it approximates a 
Gaussian in the spin labels basis, and thus satisfies the semiclassical requirements described above and leads to the 
same leading order of the 2-point function. Enhancing the ansatz for the boundary state is our second result, and 
we expect this to have implications especially in the 4d case. Furthermore, we show that it is mainly the phase term 
that is responsible for the inverse power law of the leading order of the 2-point function: the quadratic part of the 
Gaussian serves mainly to damp the remaining oscillations of the PR kernel. 

Another crucial point is the choice of the model; that is, the choice of face and edge amplitudes. An important 
question is the stability of results among different models. In fact, if one has a clear procedure to reproduce low-energy 
physics from spinfoams, then that is a powerful tool to discriminate between physically interesting models. To see 
how this can be done in practice, we consider a particular class of models, which differ from the original PR model 
only in the weights associated with the faces. Though restrictive, this generalisation is an interesting example, as 4d 
models as well have some freedom in choosing the edge weights. Here we show that changing the face weights does 
not change the leading order of the 2-point function, but only the NLO. Furthermore, the NLO is changed in a very 
simple way, and it is easy to show that whereas the original PR model (with topological face weight) minimises the 
real part of the relative NLO correction, it is the trivial face weight that minimises the absolute value of the NLO. 
This is our third result. 

Finally, we also consider a geometric issue. Both in [l] and Q , the boundary geometry is represented by equilateral 
configurations, to simplify the geometric analysis. That analysis becomes more involved for arbitrary configurations, 
but the quantities to be computed are defined in the same way, and we expect the same results to hold. To show that 
this hope is plausible, we consider a boundary geometry of a regular — though not equilateral — tetrahedron. 

The paper is organised as follows. In the next section, we introduce the toy model and describe the physical setting. 
In section ITTT1 we study the perturbative expansion of the spinfoam 2-point function, and compute the leading order 
and NLO. In section ITvl we compute the NNLO, and show that the spinfoam quantum geometry enters at this order. 
In section Ivl we discuss the boundary state, and propose a new ansatz. In section fvTI we discuss the stability of the 
results obtained against different face weights in the spinfoam model. In section IYlII we extend the construction to 
the case of an isosceles tetrahedron. 

Recall that in 3d the Newton constant G has inverse mass dimensions (in units c = 1); we define the 3d Planck 
length as £p = 16irhG, and choose units lp = 1. 



The Ponzano-Regge spinfoam model allows us to compute the physical scalar product between states of 2d geometry 
induced by 3d quantum gravity. This scalar product implements the projection on the states satisfying the Hamiltonian 
constraint, 7i \tp) = 0. The Hamiltonian constraint generates time translations, and such states define the physical 
Hilbcrt space of the theory. The standard setting is to consider a 3d manifold M representing the evolution in time 
of a 2d slice S of constant topology, M. = X x [0, 1]. We choose two triangulations di,df of the 2d slice £ on the 
initial and final boundary and a triangulation A of M which interpolates between di and df. In the following, we 
will assume for simplicity the initial and final triangulations to be identical, di = dj = d. A spin network state on 
E, \jeed), is the assignment of the SU(2) representation labels j e G N/2 to the edges e of the triangulation d. The 
kinematical scalar product is naturally constructed such that two spin network states are orthogonal if their j e 's differ. 
On the other hand, the physical scalar product is defined as: 



where the sum is over all possible assignments of irreps to the interior edges of the 3d triangulation A, and the Ponzano- 
Regge amplitude -4a (je) is defined as the product of Wigner's {6j} symbols, each associated to a tetrahedron r of 
A: 



II. PHYSICAL SETTING 
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This state sum is topologically invariant and does not depend on the chosen triangulation A. In fact, the amplitude 
defined in this way usually gives an infinite result and requires suitable gauge fixing 

The physical observables of the theory are (gauge-invariant) operators commuting with the Hamiltonian constraint, 
such as holonomy operators. These observables are constants of motion, and we refer to them as complete observables, 
in the terminology of [T^j. Using the scalar product QJ, one can compute the expectation value of a physical observable 

O, 

In the present work, we are interested in computing correlations between initial and final states for the 2d metric. 
This corresponds in a sense to computing a graviton propagator. Of course, the 2d metric is not gauge-invariant: 
these correlations will be gauge-dependent and thus a priori unphysical. Nevertheless, we expect that the insertion 
of particles in order to localize space points will allow to turn the 2d metric into a complete observable and these 
correlations into physical quantities. 

The two ingredients of our calculations are (i) working with partial observables measuring the metric fluctuations 
and (ii) using a (partial) gauge fixing defining the proper time of evolution between the initial and final slices. 
More precisely, the chosen partial observables are geometric measurements. The simplest case to consider are the 
fluctuations of the edge lengths (or of the angles) in the boundary triangulations: the relevant observable is then 
simply the spin label j e „ of a given edge eo, or its deviation Sj eo = j eo — jo with respect to some fixed reference length 
jo- As we have stressed above, this correlation is gauge-dependent. To make it a physical observable, two steps are 
required. First, we need to identify the two end points of the edge under consideration with physical points, thus 
making the function a partial observable |l5l. This is achieved by coupling gravity to a matter field: we insert a 
particles at each of these two points (see |ly ) and express the boundary state ip in terms of the field propagator (see 
SectionlVl below). The gauge fixing amounts to fixing the representation labels of some edges in the bulk A\ (9, U3/ ) 
in order to fix the (proper) time between di and df to some value T. We choose a set T of edges in the bulk and set 
(by hand) the spins j e for edges e € T to some fixed value J e {T) depending on the value of the time T. The resulting 
PR amplitude can be written as 

Al(je) = A A (je) l[S je , M T). (2) 

This obviously breaks the topological invariance of the spin foam amplitude. Nevertheless, as shown in |l4j . if T is a 
tree T (i.e. does not contain any loop), then this amounts to a gauge fixing of the diffeomorphism invariance, which 
is at the origin of both the topological invariance and the original divergence of the spinfoam amplitude. We choose 
a 1-skclcton path along A from d% to df and fix the j e 's along that path. That would correspond to fixing the proper 
time of a point particle travelling along that path. More generally we could choose a tree T in the bulk and fix the 
je's along that tree thinking once more of particles travelling along each of these edges. Finally, in the generic case 
that £ is open and has a boundary, we could insert particles on the boundary 8T, x [0, 1] and fix their proper time 
and see how this propagates to the bulk. 

Concretely, we choose a state \ip) and an observable O with which we define the correlation: 

W>\O(0)d(T)\1>) = ^(OVOOeGaJ (0V)&e8,MXOe), with Z = £ ^(j e&9i ) ^(j eg8/ ) A\ (j e ). (3) 

Ue} {je} 

As discussed in \ip) has to be a physical state of the boundary geometry, namely a solution of the full dynamics. 
Furthermore, \tp) might be a coherent state, namely a state peaked around a configuration which is a solution of the 
classical equations of motion. In the case of a free theory, for instance, it could be some kind of Gaussian state. 

We expect the two steps described above to turn our correlation into a com plet e observable. This method of 
interpreting the boundary as matter insertion will be investigated in future work |l7j : in particular, we will see how 
it determines the 3d triangulation. 

The hope is to recover the classical values of these correlations in some limit. This would have two important 
consequences. On the one hand, it would prove that we have the right semiclassical behaviour, i.e that we recover 
gravity in a long distance regime. On the other hand, it would allow us to bridge between the discrete gauge-fixing 
that we use at the spinfoam level and the gauge used to compute the correlations in the continuum case (typically the 
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temporal gauge plus the Coulomb gauge, see pj)- Clarifying these issues will definitely be very useful when tackling 
the problem in 4d spinfoam gravity. 

In this paper, we will apply this framework to the simplest case, d is made of a single (open) edge and the 3d 
triangulation A is a single tetrahedron interpolating between the initial edge and the final edge. Let us briefly recall 
the setting of Q (see also 0])- We consider a single tetrahedron embedded in flat 3d Euclidean spacetime as drawn 
in FigQ] Using the PR model, the edge lengths are given by l e = j e + \ = C(je), where we take C 2 (j) = j(j + 1) + \ 
as the Casimir of SU(2). We assume to have measured the time T = tf — 1%, and we are interested in computing the 



FIG. 1: The dynamical tetrahedron as evolution between two hyperplanes, useful to study the correlation between edge lengths. 
The labels give the physical lengths as ii = C(ji), £2 = C(j 2 ), it = C(jt), and T = tf-U = £t/V2. 

correlation between fluctuations, around a given background, of the length of the bottom edge and the length of the 
top edge. This amounts to computing a single component of the graviton propagator, such as Wn22(T), choosing 
ji and J2 to be along respectively the x 1 and x 2 axis. The measured-time setting is obtained fixing the four bulk 
edges to some representation label jt, and it realizes a temporal gauge-fixing The background around which we 
study the fluctuations is introduced by the state Assuming that \ip) peaks ji and J2 around a given value jo, we 
can compute the classical value T in terms of jo and jt ■ The generic case is described in section IVIII In the main 
body of the paper, we will restrict our analysis to the equilateral tetrahedron for simplicity and take jt = jo- In that 
case, all edge lengths are £0 = C(jo), and elementary geometry tells us that all dihedral angles have the same value 
i9 = arccos(— |), and that T = £q/\/2. 

The next case would be to consider correlations between the areas of two triangles. The 3d triangulation would 
then be a prism made of 3 tetrahedra, as in Fig|2] The temporal gauge-fixing is performed fixing the three "vertical" 



FIG. 2: The prism interpolating between two hyperplanes, and its decomposition into 3 tetrahedra. This setting can be used 
to study the correlation between triangle areas. 

bulk edges, say to a value j t ; this would leave three "diagonal" edges with unspecified labels in the bulk, which we 
will have to sum over to compute the correlations. The state \ip) will introduce a background value for the j's in the 
bottom and top triangles, say two equilateral triangles with all edge lengths £ . Wc now have nine sums to perform, 
thus this setting is computationally much more involved than the single tetrahedron setting, and for this reason we 
leave it for future investigation. We point out nevertheless that this would shed light on the following issues: 

• By allowing to compute more correlations, this setting will give the full tensorial structure of the graviton 
propagator. We could then compare it to the classical result and check the precise correspondence between the 
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discrete gauge fixing and the gauge fixing in the continuum. 

• It is the simplest setting in which we can study the effects of many tetrahedra in the bulk. 

• We would be able to analyze how the bulk edges which we have not fixed to j t get peaked (or not) around their 
classical values. This would show explicitly how fixing the time on the boundary propagates into the bulk. 

From now on we focus on the computation of the graviton propagator in the framework of the single tetrahedron, 
studying the behaviour of its leading order and developing the necessary tools to extract the quantum gravity correc- 
tions. 



III. THE FREE PROPAGATOR AND THE NEXT TO LEADING ORDER CORRECTION 



A. Computation of the leading order 



We consider the single tetrahedron model described in the previous Section. To study the perturbativc expansion 
of the 2-point function in a Coulomb-like gauge-fixing, our starting point is the following quantity, 
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(4) 



The normalisation factor J\f is given by the same sum appearing above, without the two Casimir insertions [C 2 (ji) — 
C 2 (jo)]. The latter represent the expectation values of the metric perturbation h^ v on the boundary spin network. 
\&o[ji] is the boundary state, representing the vacuum state of the theory (see the discussion in pj), and the {6j} 
symbol is the PR kernel. For the motivations leading to the expression above to be interpreted as the 2-point function, 
see 0] ■ As for the boundary state, we assume for the computation of the free propagator the following Gaussian ansatz: 



*„[7i] = exp l--6jf + iti(ji + -) j , (5) 

where Sji = ji — jo, and a = ^ (see 0)- The essential point is the 1/jo dependence of a which leads to a Gaussian 

state with a width growing in \fjo~. Remarkably, this matches the scaling of the uncertainty in measuring lengths in 
3d quantum gravity, Sj ~ y/J or equivalently 51 ~ \/l£p The exact coefficient 4/3 leads to a precise physical 

interpretation of the 2-point function as an harmonic oscillator, but changing it would not modify the asymptotic 
behaviour of the 2-point function. We will discuss later in section^the physical meaning of this ansatz. 
The large scale regime of (@J can be studied using the well known asymptotics of the {6j} symbol fToL f20| . 



Jl Jo Jo 



(S R [j e ] + f) 3 

+ 0(3 2 ), (6) 



32 Jo Jo J y/12nV(ji,j2,jo) 
where V{j\,j2,jo) is the volume of the tetrahedron, 

V{jx,h,j ) = ^/4I^H^, (7) 

with ii = ji + |, and it reduces to Vq = in the equilateral case. S^,[j e ] is the Regge action, 

fetfe] + (8) 



where the 6 e are dihedral angles, defined as the angles between the external normals to the triangles. In the Appendix, 
we give their expression in terms of edge lengths. The Regge action is a discretised version of GR, which captures 
the non-linearity of the theory. We can then expand it around flat spacctime, precisely as in the continuum. In our 
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calculation, this is induced by the Gaussian © entering ©; the expansion around the background value jo f° r both 
ji's reads 



i,k=l 



9ji9j fc i c =j 6 4 fill d 3idjhdji 



where 5r[jo] = 6$ (jo + §)■ We write this as 

S R [je] = E U^*JT«2. ( 10 ) 



a!6! 

0,6=0 



with M a b '■= a aab ■ This matrix can be evaluated from elementary geometry, and the elements which are 

relevant to the following computations are reported in the Appendix. Indeed, we have M a b ~ jo~ a ~ , thus in the 
jo 1— > 00 limit higher derivatives of the Rcgge action become negligible, and this justifies the expansion. 

As noted in the key point here is that the zeroth and first order terms in © reproduce + ]2 + 1), which 
is the phase of © induced by ©. Therefore, when we use the expansions © and © in the exact expression ©, 
the phase coming from © is cancelled or doubled, depending on the sign of the two exponentials of the cosine. But 
because the phase makes the argument of the sum rapidly oscillating, we expect only the exponential where the phase 
is cancelled to contribute to the sum. This intuition is indeed confirmed by numerical simulations, and we can thus 
write 

^1122 .70) = — -77 > — 7= . , , . exp{-iS R }, (11) 

where we have absorbed the constant factor e , 4 in the normalisation. 

4v / 3tt 

Let us comment this first result. The exact expression for the 2-point function is given in ©. By keeping only 
the first order in the large j expansion of the {6j} symbol, equation ©, we are able to rewrite © as (|llf) . This 
expression is what we would write down to compute the 2-point function using not the PR kernel, but directly the 
path integral for Regge calculus, choosing the factor 

M&) ; = np+i) (12) 

for the measure of the path integral. Notice that this measure is largely determined by the volume. The Regge 
path integral corresponds to a discretised version of quantum gravity, and it has a highly non-linear structure. If we 
now expand for large jo the measure and the Regge action, this corresponds to a discrete analogue of the continuum 
perturbative expansion of quantum gravity. The first non trivial contributions, corresponding to the free propagator, 
come from the quadratic term in the expansion © of the action, and from the trivial background term /i(jo) in the 
measure. Using the explicit value of one indeed finds Q the following leading order, 

WfSatfo) = ~ = -U\ - iV2). (13) 
2 Jo Jo 2 

This result can be interpreted as the contribution from a single tetrahedron to the free propagator ^ZT , with u> ~ 
(see Q), and it is supported by numerical calculations, reported in Fig|3|below. 

The scheme of the perturbative expansion leading to this free propagator can be summarised in the following box: 



© © {□} © D31+... 

quantum gravity — ► Regge path integral — > linearised theory 



From the scheme above, it should be clear that we have two sources of corrections to (|13|l : 

(i) contributions coming from higher orders in the expansion of the Regge action and the measure; 
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(ii) contributions coming from higher orders in the expansion of the {6j} symbol. 

We will keep this notation (i) and (ii) throughout the rest of the paper. Corrections (i) come from the non-linear 
structure of the Regge path integral; whereas the background measure and the quadratic action provide the 1/jo 
behaviour, non trivial contributions from the measure and higher orders in the expansion of the Regge action 
provide corrections in higher powers of £p. Because of their origin, these corrections can be considered as the discrete 
analogues of the QFT self-energies. 

On the other hand, (ii) are new kind of corrections, arising from deviations of the {6j} symbol from the asymptotics 
©; these corrections have no QFT analogues, and cannot be guessed by discretising GR a la Regge. In other words, 
the microscopic discrete geometry, and its non-Regge like dynamics, does affect the large scale behaviour in a non- 
trivial way. In particular, we show below that (ii) does not contribute to the NLO, but arises only at the NNLO. 

An important remark should be added here. The boundary state ^o[j] in (0J is in principle the vacuum state of the 
full theory, just as the {6j} is the kernel of the full theory. In computing the free propagator above, we assumed the 
Gaussian form JSJ. However, when computing the corrections to the free propagator, we expect contributions from 
^o[j] beyond JSJ). At the present state of investigation, the exact form of the boundary state is not fully understood, 
thus in the following we discard these terms. 1 See also the comments at the end of Section Fvl 



B. NLO: contributions from the Regge action 

Let us begin by studying the corrections of type (i); namely, we start from the Regge path integral l|ll[) and study 
the corrections arising from the measure and terms higher than quadratic in @. To make the following calculations 
clearer, let us rewrite 111 If collecting the quadratic terms in the exponent, 



31 ,32 

where we have introduced the shorthand notation 



WimO'o) = js jf X) F(jo, Sji,Sj 2 ) cxp j-i A lk 5j l Sj k ^ , 



F(j ,Sj 1 ,Sj 2 ) := Mie) II P 2 ^) ~ C2 C?'°)] ex P \ ~ l E J&i M « b 5 & 5 & \ > 

i [ a+b>3 ' ' J 



(14) 



(15) 



and we have absorbed a constant factor exp{— iAfi (jo + i)} in the normalisation. The normalisation is obtained by 
taking the same sum in (|14f) and substituting F(jo, Sji, 6j 2 ) with 



FjvUo, Sji,6j 2 ) := fi(je) cxp -i ^ ^ M ab 5j% 8j h 2 \ 

y a+b>3 ' ' J 



(16) 



Using the explicit values of a given above and of the second derivatives of the Regge action reported in the Appendix, 
the quadratic exponent can be shown to be 



Aik = a5 lk + i— 



d 2 S R 



djidjk 



1 + % cot ■& 



3,-30 3 j \ _i 1 + icQt tf 

x sin i9 



(17) 



In the following calculations, we will need the inverse of this matrix, which is 

1 e w 



ii) 



1 



Of course, we expect these contributions to be non-negligible, and the reason we do not consider them here lies in the lack of control 
we have over the vacuum state of the theory, together with the exploratory spirit of this paper. 
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+ «v2- All of our analytic calculations will keep track of this simple algebraic form. 



where e 

The fact that the dominant term in the exponent is quadratic suggests that we can evaluate the sums <|14[) by 
approximating them with Gaussian integrals. In fact, in the large jo regime, we have 52jf=o ~ Io^° dji = J"^° o d8ji 
As is well known, 



ax e 2j o 



— X 

ax e 2 ™ 



(19) 



thus approximating the sums with Gaussian integrals does not affect the NLO of <|14[) . At this point, the logic to 
study the NLO order in j is clear: we expand (|15|l and (|16|) in powers of 5j\ and 5j 2 , and then we evaluate the 
integrals as momenta of Gaussians. If we expand l|15|l. we get 



V 



(2.70 + i^SnSn + (2j + i)* Wo + n (SjiSj 2 + Sj2sh) _ 

4(jo + 1) 



, (84j 2 + 124jo + 57) ( 



n ^2 Jjjjo + 236j + 125 2 2 
+ (2.?o + l) on , . , 7Tn (ojidJa +dj 1 dj 2 ) + (2j + 1) 1 . , 0j 1 dj 2 



32(j + l) s 



16(70 + l) s 



a+b=3 \a+b=3 



8j\ 



(20) 



The structure of this expansion is rather intricate, because the coefficients go down as jo <—> oo, whereas the integrals 
(approximating the sums) go up. In more detail, we have that on the one hand, the coefficients of Sjf in the expansion 
of the measure go as — i=j 3 ~ i: , and M a b ~ jo~ a_b ; on t nc other hand, because the squared width of the Gaussian, 



and thus the values of the squared momenta, are proportional to jo , each factor Sjf in the integrand contributes a 
factor j after the integration is performed. By inspecting (|2U)l . we see that the leading term is of the form jo/\/Vo, 
and it comes from the first term in both the expansion of the measure and the integral. The next to leading is of the 
form jg/VVo, and it has a number of contributions: terms up to the fourth order in the expansion of the measure, 
together with the first, second (with a + b = 3, 4) and third (with a + b = 3) terms in the expansion of the exponential 
contribute. Namely, terms up to the fourth order in the expansion of both the measure and the Regge action enter 
the NLO. 

The same analysis for the expansion of the normalisation J\f shows that we need only the second order in the 
expansion of the measure, and still the fourth in the expansion of the action, 



-FV(io,<5ji,<5j2) 



1 



V 



(2.70 + I) 2 - i(2j + l)(Sh + Sj 2 ) + ^(Sjf + tfjf ) + ^jitfa 



1 ~ i E jfci M - *& - \{ E jfr? M <* *if 

a+b=3,4 \a+b=3 / 



(21) 



At this point, we evaluate separately the Gaussian integrals in the numerator and in the denominator. We then 
collect all terms in decreasing powers of jo, thus obtaining the following structure, 



Wii22(jo) 



1 hjl + b 2:] l + ... = J_ 
j 4 /;,./,' + B 2 jo + ■■■ n 



— + —(b2 - 
Jo Jo B i 



(22) 



3e l 



In particular, h - 
be now computed to be 



B 1 = 



and gives the LO <|13[) discussed in the previous Section. The NLO can 



Vdet A 



TX/-NLO (i) 

"1122 



1 7481 - il048V2 1 , 

il 5184^ — -2(1-4431-0.2859), 



(23) 



where we have again used -d = arccos(— 1). The main point to remember is that this is the contribution to the NLO 
coming from corrections of type (i). 
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To support this analytic calculation, wc fit the NLO using the numerical evaluation of the exact formula and 
subtracting the fit for the LO (|13fl . The result of the fit is a behaviour with numerical coefficient —1.4431 + i0.2860 
(see Fig|3J|. This matches our analytic calculation up to three digits. Therefore, there is room for corrections of type 
(ii) at NLO only if these are very small, of magnitude 10~ 4 . 




FIG. 3: The plots showing the numerical evaluation of the 2-point function. Left plot: the leading order. In red, the numerical 
evaluations of the dashed green line is the analytic result 1131 1. Right plot: the leading order is plotted again, this time 
together with the next to leading order. The numerical evaluations of the latter is in blue, whereas the analytic result i'2'M) is 
the red dashed line. 



C. NLO: contributions from the {6j} 

We now consider (ii), the contribution coming from higher orders in the asymptotics of the {6j} symbol ©. From 
the considerations at the end of the previous Section, we know that if these corrections are present, they are indeed 
very small. We begin the investigation by studying numerically the next term in the asymptotics of the {6j} symbol 
for the equilateral case, when all spins are jo and SrIj'o] = 6(jo + We use the following ansatz, 

and we perform a two parameters fit, finding the values ai = —0.3575, ip = 0.6843. The good performance of this fit 
suggests the simple ansatz 

q 2 cos(S R [7i, j 2 , jo] + <p) , . 

(6^2)1 V(h,j 2 ,jo)i 

as the second term in the expansion ©. If we assume this ansatz, we see that the second term in the asymptotics 

goes as j . This is just one power below the first term, thus in principle we should expect the second term of the 
asymptotics to contribute to the NLO of Wn22(jo)i which is precisely one power below the LO. Notice also that a?, is 
of the same magnitude as 2170^ ~ 0.4744, so that if the second term does contribute to the NLO, it will do so with 
the same magnitude as the contribution (i). However, we know that there is no room for this, as the value computed 
from (i) already matches the exact numerical value up to 10 -4 . 

We conclude that corrections coming from the second term do not enter the NLO, i.e. 

= W$g l \ (26) 
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Understanding the exact way this happens would require a more complete study of the second term in (|24[1 . which is 
beyond the scope of this paper. 2 



IV. NNLO: THE PURE QUANTUM GRAVITY EFFECTS 

To study the NNLO, we proceed as in the previous section. We first compute analytically the corrections of type (i) 
coming from expanding the measure and the action, and we then compare the result with the numerical fit obtained 
from |@} , to see if corrections of type (ii) enter. The analytic calculations can be performed as above, thus we do 
not report here the somewhat tedious details. However, the complexity of the calculations grows as the number of 
contributions increases largely with every order; we notice that for the NNLO, in the expansion (|20f) . we need up to 
the eighth derivative in the expansion of the measure, and up to the sixth derivative in the Rcgge action (namely 
a + b < 6); similarly, we need up to the sixth derivative in both the measure and the action in the expansion of Fj\f . 
The result we obtain is 



= \ 62758 t;Z 789V ~ 2 - ^(0-2522 + ,0.7263). (29) 
il 248832 $ K 



As expected, the NNLO is one power below the NLO. Notice that the magnitude of the real part is one order of 
magnitude smaller than the NLO, and the magnitude of the imaginary part same order. This is going to be changed 
by (ii), the corrections coming from the quantum gravity kernel. 

As before, we now evaluate numerically the NNLO, by subtracting H13|) and (|23|l to Q. The best fit gives 

W NNLO( t o t) = ^NNLo,;, + j^ofl „ J_ (_ . 02 25 + iO.0730) . (30) 

Jo 

This fit is quite different from (|29|1 . which means that the corrections (ii) do enter the NNLO. Notice that, as expected 
from (|24|) . (i) and (ii) are of the same magnitude. What is remarkable, is that (ii) have the neat effect of reducing the 
magnitude of the correction. At this point it becomes interesting to evaluate analytically the second order in H24JI . 
and then compute exactly the NNLO coming from it, W^ 1 ^ '"'- We leave this issue open for further development. 

Let us now use this result to draw some conclusions. Our analysis of the NLO has pointed out that it arises only 
from higher order corrections in (| 1 1 ft . Namely, the NLO is a result of the Regge path integral alone, and is not 
sensitive to deviations arising from the spinfoam quantum geometry. The latter, on the other hand, are important at 
NNLO, where they contribute with the same order of magnitude as the Regge-like corrections, but with opposite sign, 
so that the total NNLO coefficient is smaller than that of the NLO. The structure of the corrections is summarized 
as follows: 



2 To give a feeling for how this might happen, we note the following. Using the ansatz 1251 . which is the simplest homogeneous extrapolation 
of the second term in 1241 . the expansion around the equilateral background would be 

cos (Sr[3'i, j'2, jo] + f) _ a' cos (5 R [jo] + f) 



(12 



(6^2) ( 



V(jl, j 2 ,jo) 6 



30 



Vn 



(27) 



where a 1 = J a 2 . The expansion 1351 would now look like 



W1122 (jo) 



1 6lJ§ + 



b' 2 j A 



3o Biil + B 2j'o 



1 



bi 1 / / ; 

— + ^(b 2 + a'b' 1 - (B 2 + 

30 Jo v 



(28) 



It is easy to see that 1271 does not contribute at NLO, but only at NNLO. In fact, the effect of this new term simply amounts to a 
change in the measure 1121 . thus 6' x = b\ and B' ± = B\\ therefore the NLO of 1281 is ^-(^2 — sT^i)) unchanged by the addition of 
1271 . On the other hand, b' 2 ^ b 2 and B' 2 ^ B 2 , because these terms are sensible to different measures; therefore 1271 contributes to the 
NNLO. We will see in the next section that this is indeed the case. 
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full propagator 



Rcgge path integral corrections 



W = VK LO + W njj ° + W 



NNLO 



+ 



quadratic term in the 
Regge path integral 



Regge path integral corrections 
and 

spinfoam quantum geometry corrections 



It would be very interesting if a similar structure emerges in 4d as well. To this regard, we conclude this section with 
a speculative remark. It has been argued that background independent approaches could be the key to avoid the non- 
renormalisability of GR. The PR model considered here is one such approach, and 3d GR, if treated perturbatively by 
expanding g^ u = ?/ M „ + h^ v . suffers the same infinities as the 4d case (in a gauge where h^^o). Thus, it makes sense 
to address this idea within the 3d case. In a discrctised picture, the conventional pcrturbativc framework is equivalent 
to using only the Regge path integral, namely starting our computations from i|ll|) . So, in a sense, the perturbative 
expansion that includes only corrections of type (i) is the discrete equivalent of the conventional non-renormalisablc 
perturbative expansion. However, because we did not start from Ullfl . but from the exact background independent 
quantity Q), we find that new corrections will come in at higher orders, thus changing the structure of the perturbative 
expansion. Could this provide a possible mechanism for renormalisability? It is very speculative to say so, however 
we find it intriguing that the new corrections enter at NNLO, precisely where the non-renormalisability of GR shows 
up. 3 



As should be clear from the calculations above, the boundary state plays a key role for the extraction of the inverse 
power law behaviour, thus it is important to motivate the choice |J5J|. In principle, the boundary state should be given 
by the vacuum state of the full quantum theory; however, such a state is not known in LQG. To overcome this lack 
of knowledge, one can proceed with an ansatz. To compute the free propagator, in Q it was proposed to take, as the 
vacuum state for the linearised theory, a Gaussian state peaked around the values of the labels for both the intrinsic 
and the extrinsic geometry of a chosen background, by analogy with quantum mechanics. This suggestion is realized 
in 3d by JSJ). Indeed, we can think of it as a state for a quantum mechanical system of conjugate variables j and 
</>, representing discretized intrinsic and extrinsic geometry. Then, the main property satisfied by JSJ is to peak the 
values of those variable around, respectively, jo and and to minimize, in the limit jo <— > oo, the relative uncertainties 
and ^fj^ . Notice that in order for this last property to be true, it is crucial the 1/jo dependence of the inverse 



Below we show how this construction can be made more precise, by looking at j and <p as conjugate variables 
as defined by the Fourier transform on the structure group of the PR model. Before doing that, we discuss some 
interesting aspects of the phase term in (JSJ). 



V. 



ON THE BOUNDARY STATE 



width a. 



A. The phase term 



We start our investigation by noting an important fact: it is the phase term of JSJ 



(31) 



3 In |21| it has been shown that the spinfoam quantum geometry indeed affects the propagation of point particles. Interestingly, it does 
so in such a way that the new dynamics can be encoded in new Feynman rules. It would be interesting to study if the modifications 
produced here can be as well reabsorbed in modified Feynman rules. 
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that is mainly responsible for the 1/ jo behaviour of the asymptotics. To understand this, recall that the phase term of 
represents a state whose dihedral angles are peaked around the equilateral configuration. Therefore, the presence 
of the phase term is enough to expand the Regge action around such a configuration; but since the bulk edges of the 
tetrahedron arc fixed to jo, the equilateral configuration is exactly the one used in the calculations in Section lnTl This 
means that even if we set a = in (JSj, we can proceed consistently with the perturbative expansion as in Section 
IIIII The only difference is that now the matrix A in l(T7|l consists only of the second derivatives of the Regge action. 
Indeed, we find the following analytic result: 



1 9 

W / 1122(jo) - — i 



1 



30 2V2 3o 



245 
~64~ 



■ 9 ^ 



(32) 



Once again, we support this result with the numerical simulations, which are reported in Fig^J Notice that we get 
much more oscillation, due to the absence of the quadratic damping factor. In particular, we see from the right figure 





FIG. 4: The plots showing the numerical evaluation of the 2-point function with a — 0. Left plot: the leading order. In red, 
the numerical evaluations of @; the dashed green line is the analytic result in 1321 . Right plot: the numerical evaluations of 
the next to leading order, in blue, together with the first two terms of the analytic result 12311 . The green dashed line is the 
leading order, whereas the red dashed line shows the next to leading. 

that the NLO oscillates quite strongly, and this makes it hard to fit its value. Of course, JSJ with a = loses the 
semiclassical properties described above, however, this is just an exercise to show that actually is the phase term the 
key object, as much as the calculations are concerned. 

The fact that the phase term alone is sufficient to get the right leading-order asymptotic behaviour, suggests that 
we have freedom in choosing the angle: with a = 0, different angles would still give non trivial results. Indeed, these 
would correspond to non equilateral geometric configurations (see Section IVIII below 1 ) . 

The phase term has another interesting property. By inspecting Q, we see that it is the only (D-valued quantity 
entering the definition of the 2-point function; thus, it is responsible for the 2-point function being a complex quantity. 
Recall that we are working in a measured-time setting, where the 2-point function is reconstructed 
complex result is what we need in order to be consistent with this interpretation. However, one can also consider a 
different setting, where all boundary edges are allowed to fluctuate. This "general boundary" context is particularly 
relevant in the 4d case model discussed in Q], where the measured-time setting is not applicable. 4 In the general 
boundary context, one expects to reconstruct the 2-point function in configuration space; in the Riemannian case, a 
real function of the spacetime distance between the points, so that the computation of the 2-point function has to 



4 This can be easily understood, as the setting described in section [n] (see in particular Fig0 does not generalize to 4d: there is no 
embedding of a 4-simplcx in 4d Euclidean space such that two triangles of the 4-simplex lie in different parallel hyperplanes. 
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lead to a real result. This is not obvious, if the phase term is (D-valued. In fact, notice that in 0, the author has 
to choose a (D-valued (non-diagonal) quadratic term in the Gaussian boundary state, in order to compensate for the 
complex phase term, and make the final result real. Here, we would like to suggest that the reality of the result can 
be obtained keeping the simple quadratic Gaussian term, but rather changing the phase to a real quantity. We argue 
below that this can be done without losing the 1 /jo behaviour. 

The immediate choice to make the phase real is to take either the real or the imaginary part of <|31(l . i.e. e'^+s) 15 ± 
e -«(j+2) 1 '. In principle, it is not obvious that such a phase term will still guarantee the 1/jo behaviour, thus we find 
it useful to address this issue within the toy model so far described. To fix ideas, let us choose the real part, and 
consider a new boundary state in Q whose phase term is given by Yii cos(ji + ^)$. We keep the same quadratic term 
as in and we use again © to study the leading order. This leads to a structure of the type 

Recall that the linear order of the Regge action is [4(jo + \) + ji + jz + l]i9. As discussed in Section ITTT1 only the 
exponential of the Regge action when the linear term is cancelled will contribute to the 2-point function, the other 
term being negligible due to the oscillations in the sum. It is thus clear that taking the cosine as the phase term 
makes in principle both exponentials of the Regge action matter. However, attention must be paid to the factor 
4(jo + h)^ in the linear order of the Regge action. In Section ITTT1 we were able to discard the related phase e -44 ^ ^)' 5 , 
by reabsorbing it in the normalisation. This can not be done now, because both orientations matter. Including also 
the normalisation, wc indeed have 



^1122 (jo) 



J 4 e - l4 U + i)^- 4 f [ Blj 2 + ...]+ e i4(jo + *)i»+iS [ Blj 2 

1 Rc [b l3 l + ...]+ Jm\big + ...] tan [4(j + |)g + if] 
Jo 4 Rc [B^l + ...]+ Im{B l3 % + ...} tan [4(j'o + |)0 + *5 ] ' 



(34) 



While the first terms of numerator and denominator would produce the right 1 /jo behaviour, we see that we have an 
extra jo-dependent phase which spoils the result. Therefore, the choice Yii cos(ji + ^)"& does not work. It is easy to 
see that this problem can be fixed by taking 

= cos(j ?; + 2j + (35) 
so that the linear term of the Regge action is exactly cancelled. In this way we get 

W c\ 1 Re [6 ie -'f jQ 5 + . . .] 

W1122U) = — p r D _ ilL , 2 : r (36) 

Indeed, this latter choice cancels all jo-dependent phases, and the 1/jo behaviour is restored. In particular, using the 
same quadratic term as in JSJ), we obtain the same b\ and B\ of (|22|l . and thus the LO of (|36|l is 

__, ... 1 / Rebi +Im6i \ 3sin|z? 1 

This result is indeed supported by the numerics (see FigEJ. Therefore, it is possible to keep the simple diagonal and 
real quadratic term in the boundary state, and choose a phase term such that the 2-point function is a real quantity 
at all orders, with leading order proportional to 1/jo- 

However, let us stress that this construction is not physically interesting in the measured-time setting, where the 
2-point function is a complex function, as mentioned above. This is reflected in the fact that (|35|) is not a sensible 
choice for the phase term: the boundary is defined on the hypcrplanes and it should not refer to the bulk edges. 
On the other hand, the interest in the computation performed above is to have an exploratory investigation towards 
the general boundary context. In fact, in the general boundary setting, all edges e of the tetrahedron are allowed to 
fluctuate, thus it makes sense to take rie cos (j e + as the phase term of the boundary state. Since this quantity 
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FIG. 5: The plots for the real boundary state l|35|l . Left plot: the leading order. In red, the numerical evaluations; the dashed 
green line is the analytic result in 1371 . Right plot: the numerical evaluations of the next to leading order, in blue, together 
with the analytic result of the next to leading in 1371 . 



contains a phase that cancels exactly the linear order of the Regge action, we expect its use to give a real 
behaviour. The calculation performed here supports the idea that such a choice would give a correct real 2-point 
function. In the measured-time setting, the only interest in (|35J) is to show that using a real phase term would still 
work, with the benefit of a real 2-point function. 

We conclude with a remark about choosing the imaginary part of (|31|l . Because of the normalisation in Q), this is 
equivalent to take the SU(2) character itself, i.e. 

<Mi) = tt(o) = ■ (38) 

2 2i sin j 

It is interesting to note that (|38J) can be interpreted as the propagator of a point particle coupled to the PR model [2lJ . 
Consequently, the use of (|38[) as the phase term of the boundary state opens the way to a new possible interpretation 
of the gauge-fixing performed in this computation: (|37[) can be thought of as the 2-point function in the gauge fixed 
by the presence of a particle in the boundary. 5 



B. Harmonic analysis 



We parametrise SU(2) elements as 

g((f), n) = cos ^ 1 + i sin ^ h ■ a, <f> G [0, 47r[, h € S 2 , 

where oi are the Pauli matrices, satisfying of = 1. Moreover, the group element <?(</>, n) is obviously identified to 
g{—(f), —n). We can thus restrict 4> to live in [0, 2tt]. Finally, g((f>,h) and g(4> + 27r, n) = —g(<p,n) actually define 
the same rotation in the 3d space (same SO(3) group element). In this representation, the characters and the Haar 



5 With this interpretation, it is natural to suggest that the dihedral angle 1? is the deficit angle produced by the particle on the boundary. 
But recall that the explicit value of the dihedral angles of the boundary depends on the choice of triangulation, •& being the case of a single 
equilateral tetrahedron; thus, this interpretation suggests a possible link between boundary matter and finiteness of the triangulation 
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2 ■ 2 ip 

FIG. 6: We plot the three considered Gaussian distributions on SU(2) peaked around ±1: fi(4>) = e~~ sln 2 , f2{<j>) = 
2~2j Xj(0) e ~" 2 ""' an d f:j(<t>) = e - ^ 1 * (for £ [0, 7r] and its mirror image e - ^^"" - **' for (f> £ [7r, 27r]). The two plots respectively 
correspond to the values a = 2 and a = 0.2. On the first plot, it is clear that the three Gaussians are distinct and that /3 is 
not smooth. On the second, we see that /2 converges toward the distribution /3 as a goes to (high jo), thus justifying 139t 
while the distribution f\ dual to the Bessel functions remains slightly different. This deviation of /i from the standard /2 will 
be responsible for the different NLO of the 2-point function when using it as boundary state (see below). 



measure are given respectively by: 



siiuj + ±)0 _ 1 



= sinf ' dg = 4^ Sin 2 d ^' 

A class function / is a function on the group which is invariant under conjugate action, f(g) ~ /(hgh" 1 ). With the 
paramctrisation chosen above, it is a function of the angle cf> only. For such a function, the Fourier transform can be 
simply written as /(<?) = J2jfjXj(g)> with fj = J dg f(g) Xjid)- If j an d 4> in © are conjugate variable like this, 
then just as (JSJ is a Gaussian peaked around jo in the irrep labels basis, then its Fourier transform should describe 
a Gaussian on the group peaked around the class element However, in spite of its simple structure, © does not 
admit a simple Fourier transform. For, in the limit jo > oo, we can approximate the sum with integrals, and find: 

*oW = E*o(j)xi(0)=E e " f(J '" io)2+iO ' +4) "^(0) 

3 

_t L*(io+i)(tf+<« e -2t(0+tf) 2 _ ^Oo+^Xtf-^g-^W-tf) 2 ! _ (-39) 

sin | L 

However, notice that a function of the type e - ^^ is not smooth (or even continuous) on SU(2). This is due to the 
compactness of SU(2) and the gluing condition at the "boundary" = [2ir], 
A smooth, more natural Gaussian-like function on SU(2) would be 

Mg) = = j^e-^ 2 K (40) 

where the trace Tr is taken in the fundamental j = 1/2 irrep. This function can be easily decomposed into its Fourier 
modes; notice in fact that Tr(<7 • a) 2 = x|(<?) — 4, thus (|4(J|> can be expanded in a simple sum of characters, thanks 

2 

to the decomposition property of tensor products into irreps, xl = c jXj- This feature makes it a more natural 

2 3 

object to use within the group structure of SU(2). For a comparison with (|39(l . see FigH3 
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To find the normalisation N' , we evaluate the integral 

,-2k 



| TV' 



dg e 

SU(2) ' n 







dd> sin — e 
Y 2 



Io(-)-h(-) 

a a 



(41) 



where I n (—), n € N, are the modified Bcsscl functions of the first kind (see the Appendix). In a similar way, we can 
compute the Fourier transform of 14Ufl . For j € N, we have 



M) = 



SU(2) 



dg-ipo(g)xj(g) 



l l 



d<p sin — sin(2j + 1)— e 



sin 2 # _ ~ 



(42) 



where we used the explicit value l|41|l. On the other hand, ipo{j) vanishes if j is strictly an half-integer. 

At first sight, it is not obvious that (|42H fulfills the semiclassical properties required earlier. However, as we show 
in the Appendix, the a i— > limit of (|42|) is \[2 \f^&~^^ , namely a Gaussian peaked around j = 0. We are on the 
right track. To obtain a function peaked around jo, we need to shift the peak. This can be done in a simple way, 
now that we understand the Gaussian in the irrep label basis as the Fourier transform of l|4(J|) ; to shift the peak of the 
Gaussian to a value jo, it is enough to multiply (|4(J[I by the appropriate Fourier mode, namely Xjo(d)- Accordingly, 
we consider the following function, 



Ma) = le^<»-*> s 



Xj (g), 



N 



Mr)-W(i). 



(43) 



Proceeding as above and using Xj(g) Xjo(g) = J2k Xk(g), it is straightforward to check that the Fourier transform of 
(g3J) is 



Mj) 



1 

N 



3+ jo 

E 

k=\j—jo\ 



dge^-*)\ k (g) = T ^°^0 

SU(2) /t„(2\ r_ . . ! 2 ' 



(44) 



Again, this result holds for j integer, whereas ipoti) vanishes for j strictly an half-integer. For small a, this is 
approximated by (see the Appendix) 



V'o(j) 



J |j-3o|(^) ~t/2*[® ( 



rO'-io) 2 



(45) 



and it describes a Gaussian peaked around jo, with squared width i, as we wanted. 

Let us now discuss the dihedral angle. I|43|) is a Gaussian peaked around = 0. We want to shift this peak as 
well, around <f> = This can be done adding a phase in j. The simplest way to do so, without spoiling the property 
of having a simple Fourier transform, is to add to l)44|l the real phase term cos(j + |)$, which we discussed in the 
previous section. Then, we consider 



= ^T'^ cos(j + 



(46) 



The reason why the cosine is the easier choice to add a phase in j is that using 



cos 



3 + 



1 



i) 



xM = 



2 sin | 



2 sili I 



xj (<!> + #) 



we can simply compute the Fourier transform of Ij46(l shifting the character's argument as in l|47jl. and we get 

1 



*o(ff) 



2N sin § 



sin(^) X30 ^-^)e- 



^ + sin(^ )x , (^ + ^)e-*- 2 



(47) 



(48) 
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Let us comment on this expression. In the general boundary context, l|4bj) can be taken as a new ansatz for the 
boundary state. It has a good semiclassical behaviour, in the sense described above: in the jo >— > oo limit, it becomes 
a Gaussian peaked around j = jo (see (|46fl 'l and </> = ir — (see 114811 s ). Furthermore, the interpretation of jo an d $ 
as conjugate variables has now a clear mathematical meaning: they are conjugate with respect to the SU(2) Fourier 
transform. 

Notice that we are led to include also the configuration <f) = ir — d: the SU(2) angle <f> represents both the external 
dihedral angle •& and the internal angle 7r — $ between the triangles. In other words, the boundary spin network has 
an inside-outside symmetry. 

In the general boundary context, we expect the use of (|48|l to lead to a real propagator W with the correct 1/jo 
leading order. This is suggested by the asymptotics (|45f) . To support this expectation, we consider this new ansatz in 
the toy model presented in this paper. In our context however, we cannot use the cosine phase term, as we discussed 
above. Then, we consider (|42[) for the quadratic term and the the usual phase (ST} . Choosing the same a = ^ as 
before, we obtain from the numerics exactly the same leading order l|13[) (see also FigQ, 

W1122O0) -(-0.5 + HA) + 4(-1.3 - »0.2) (49) 
Jo Jo 

Notice that the next to leading order is now different. 6 This is due to the fact that l|44|l contributes non trivial 
corrections. Indeed, the leading order of its perturbative expansion, 1)45(1. reduces to the quadratic part of J3J|, 
but higher order terms enter the calculation of the propagator corrections. From this point of view, the use of 




FIG. 7: Leading and next to leading orders using l|44|l . Left plot: the leading order. In red, the numerical evaluation; the 
dashed green line is the analytic result in <\'A'2l . Right plot: the numerical evaluations of the next to leading order, in blue, 
together with the leading order plot. 

(|44|l enriches the toy model, and it makes it possible also to consider the most complete structure of the graviton 
perturbative expansion, where both the kernel and the boundary state are given as a power series in £p. 

As for the contribution of these new corrections, notice that the NLO in (|49|1 is slightly smaller than the NLO J22J). 
This might look like a side result. However, in 4d the numerical complexity grows very fast, a fact that has slowed 
down numerical computations thus far. Therefore, having a new boundary state which reaches the same asymptotic 
behaviour earlier could turn out to be very useful, were the same behaviour to hold in 4d. Furthermore, it turns out 
that using such a boundary state allows an exact calculation of the LO of the 2-point function in the 4d case; namely, 
one does not need to use the (troubled) asymptotic formula of the 1 10 7 1 symbol, but one can uses its exact expression 
in terms of integrals of character. A work on this is in preparation [23j. 



The reason for fewer digits of precision lies in the numerical difficulty of calculating the Bessel functions. 
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VI. COMPARING DIFFERENT SPINFOAM MODELS: STABILITY OF THE LEADING ORDER 



The definition of the n-point functions in the spinfoam formalism is not only fundamental to extract physical 
predictions, it is also important for distinguishing between the different existing models. In particular, obtaining the 
correct LO for the 2-point function is a requirement for any sensible model. In 4d, there are a variety of different 
models, and using the equivalent of to study their large scale limit could possibly be a good way to narrow them 
down. However, some variety is also present for the 3d models, and it is useful to consider the 3d toy model to begin 
to investigate how different models can produce different 2-point functions. 

We consider a simple generalisation of the PR model, where we modify the face weights, by taking 



•A k A (je)=U(2j e + l) k l[{6j} T 



(50) 



as the amplitude. The original model is reproduced for k = 1. 

From the point of view of the Regge path integral (|llfl . this new definition of the kernel amounts simply to a change 
in the measure term (|12|l . which becomes 



n (2ji + i) fc 



(51) 



If we take into account this change, we can straighforwardly recompute the LO and the NLO, and we find that the 
fc-dependence enters in a very simple way: 



1/1 1/1 



2161 
1296 



7 



:fc + 



17 



3V2 648%/2 



)]• 



(52) 



Let us comment on this result. Firstly, the LO is insensitive to this kind of modification. This means that any spinfoam 
model of the type (|50|l leads to the same asymptotic behaviour of the 2-point function. This result is expected, as we 
know from the previous analysis that at leading order only the trivial background value for the measure enters, so 
that any fc-dependence is washed away by the normalisation in Q . 

On the other hand, the NLO depends on k, thus different models lead to different corrections. While the (conven- 
tional) case k = 1 remarkably minimises the real part of the relative correction JT NLO /JT LO (see left plot of FiglHJ, 
the minimum of the absolute correction is reached for fc = 0, the case of trivial face weights. 7 This feature is clarified 




4- 




3- 












-1 


1 




k 



FIG. 8: The relative correction W NLO /W LO , as a function of the measure weight fc. On the left, the case with boundary state 
0; the quadratic function is the real part, while the linear function is the imaginary part. We see that fc = 1 is the integer 
value closest to the minimum of the real part. On the right, the real case with (1351 as the phase of the boundary state. We see 
that fc — is the integer value closest to the minimum. 



7 Interestingly, it has been advocated in the 4d case that a model with trivial face weights could be better behaved |24| . 
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if we use (13811 as the phase of the boundary state, which leads to a real 2-point funtion. In this latter case, the 
/c-dependence of the correction is 



rk fs s 1 1 /3 ; o 5 , 8701 



WTiaatto) = -— - -(-e + —k + T7 ^), (53) 



2j Q j^A 12 15552 
and we see that the case k = minimises the correction (see the right plot of Fig|SJ) . 



VII. EXTENDING THE TOY MODEL: THE ISOSCELES CASE 



In this last section, we want to address a different issue. In the toy model considered so far, we have restricted 
our attention to the case when the background is chosen to be an equilateral tetrahedron. However, we expect the 
results here described to extend to a generic background configuration. To give an idea of the way things would work 
out, and as a first generalisation, we consider in this Section the case when the two opposite edges of the background 
tetrahedron have spin jo, and the four bulk edges have spin j t . We call the corresponding dihedral angles -do, $t (all 
relevant formulas are reported in the Appendix). In this case, we can still consider an embedding in 3d Euclidean 
spacetime as in Fig^ an d T = t 2 — t\ = \J 2jf — j 2 . 

The physical setting of the problem is exactly the same. The only difference is the geometric analysis, which is now 
slightly more involved. In particular, notice that not all configurations (jo,jt) are admissable; since the tetrahedron is 
embedded in Euclidean spacetime, triangle inequalities must hold, and thus we restrict in the following to \/2jt > Jo- 

Using the simple ansatz of the type JSJ for the boundary state, the expression Q for the propagator now reads 



W 1122 (T) - 



1 1 



2jt 



\{2j i + l) [^(jO-C^jo)]^ 





Jl 


") 


\ 32 


3i 


It / 



(54) 



with 



a 1 
*o[7<] = cxp { --{ji - jo) 2 + iMk + 2) 



(55) 



and cos$o 



4j? 



(see the Appendix). To obtain the leading order of l|54|l . we study the limit (jt,jo) ^ 00 



with r = jt/jo fixed. This last requirement is necessary in order for © to hold; furthermore, it would not be very 
interesting to study the inhomogeneous limit j t <—> 00 with jo fixed, because this means going towards a degenerate 
configuration where the tetrahedron is squeezed into a 2d surface. In the limit (jt,jo) l— * o°j we proceed as in Section 
IIIII using © for the {6j} symbol, and expanding the Regge action. Choosing the value 



4r 



a 



jo (4r 2 



(56) 



for the free parameter entering l|55|l . we obtain 



W 1122 (j u j ) = ~£ F (h, 631,632) e" lE AtkSn&Jk 
3o -A/ tt 



(57) 



with 



A 



ik 



jo(4r 2 



4r 



1 + i cot ■&* 



sin i?t 




(58) 



Here we used the fact that cos^t 



(see the Appendix). Once again, this matrix can be interpreted as 



the kernel for a harmonic oscillator (see for instance Q, or the appendix of 0), if we identify the frequency of the 
oscillator with u = Indeed, it is this request that motivated our choice JS5J|. 
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At leading order, F(Jq, Sj\, 8J2) can be read from the first term of l|2(J|l : the analogous term in the normalisation is 
the first term of l|21|l . Approximating the sums as Gaussian integrals as in Section IlIII we get 



Wi 122 (jt, jo) 



4r 2 - 1 
2jt 



(59) 



For r = 1, we recover the equilateral result l|13|) . To consider non equilateral configuration, we plot H59|) for different 
values of r. In particular, if r ^ 1, we have cosi? t ~ | and w ~ and (|59"j) reads 





FIG. 9: Exact numerical evaluation of IW1122I as a function of jt, for r = 2.5,25, and 50. On the left, the leading order, 
computed in 15911 . On the right, the next to leading order, not computed analytically here. 



2-7T p iulT 

Wll 22(T)c- — . (60) 
Jo 2w 

A comparison of the exact numeric calculation with the asymptotics is shown in Figure El where we see good 
agreement with the prediction of (|59|l . Note, however, that the asymptotic value is reached more and more slowly as 
r increases. 



VIII. CONCLUSIONS AND PERSPECTIVES 



We have considered a single tetrahedron embedded in 3d Euclidean spacctimc as a toy model to test how the 
spinfoam formalism can be used to construct the perturbativc quantum theory of gravitons. We focused on a single 
component of the propagator, defined in and computed the NLO and NNLO corrections to the free l/£ behaviour. 
These are given respectively in <|23")l and (|29"j) and, as expected, have the behaviour l/l 2 and l/£ 3 . 

We found that these corrections have an interesting structure. On the one hand, the NLO can be entirely computed 
from an expression, given in l|ll|) . which is a path integral for the Regge action. This can be seen as the discrete 
analogue of the conventional continuum definition of the 2-point function. On the other hand, the NNLO gets 
contributions from the true, non Rcgge-likc spinfoam dynamics. This in turn cannot be written down by simply 
discretizing the GR action a la Reggc. Here is where the use of the spinfoam formalism modifies substantially the 
conventional expansion. It is particularly intriguing that this new source of corrections enters at NNLO, precisely 
where the continuum theory becomes non renormalisable, and that it reduces the magnitude of the NNLO correction. 

We introduced a new ansatz for the boundary state, given in (|46ll . which reduces to the original Gaussian in the 
large j limit, and thus satisfies the semiclassical requirements, but has the advantage of making j and the dihedral 
angle conjugate variables with respect to the harmonic analysis over SU(2). This property helps treating the sums 
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appearing in the evaluation of expectation values within the PR model. In particular in 4d, the Riemannian Barrett- 
Crane model uses the gauge group Spin(4) ~ SU(2) x SU(2), to which we can apply all the tools developed here 
for the gauge group SU(2). Remarkably, the use of l|46|) in the general boundary context allows us to write the 4d 
propagator as an integral over the group ■ 

It would be interesting to extend these results to the Lorentzian case, using the expression Q with group SU(1, 1). 
The SU(1,1) {6j} symbol satisfies the Lorentzian version of the asymptotics I© [2(|, thus we expect to be able to 
perform the same analysis pursued in this paper for the Riemannian case. It is also worth mentioning that in the 
non-equilateral Lorentzian case, there are no triangle inequalities to fulfill, thus the analysis is richer. 

In light of these results, we think that this approach is very promising: it allows computation of the leading order as 
well as the corrections of (a component of) the graviton propagator, and, in a more general perspective, of scmiclassical 
observables. When applying this program to 4d, it will be interesting to compare and possibly match the spinfoam 
corrections to the loop corrections computed in the standard perturbative expansion (see e.g. Q). Let us nevertheless 
point out that the NNLO in our 3d model already contains an impressive number of terms and that the technical 
difficulties to organize and compute further spinfoam corrections will likely be comparable to the ones encountered 
when computing and summing the Feynman diagrams at two-loop and higher orders. 

Important issues remain open. The next step, to our opinion, is to be able to compute all the components of the 
propagator, and recover the full tensorial structure. This is crucial in 3d, in order to check satisfactory that the theory 
does not have local degrees of freedom. The other fundamental question is the definition of a suitable coarse-graining 
procedure. This is needed to understand how the single tetrahedron contribution is related to the true large scale 
limit. Both these issues require the extension of this formalism to deal with many tetrahedra. We sketched how this 
can be done in Section II, and we leave this open for future work. However, it is important to realize that increasing 
the number of tetrahedra makes the numerical evaluations much more difficult. 
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APPENDIX A: SOME ELEMENTARY GEOMETRY 



In this appendix we report the geometric formulas used in this paper. We begin by the derivatives of the Regge 



action Sr, namely the matrix M a b 



dS R 



3e=]0 



introduced in © . Recall that i e = j e + i , and 



(Al) 



All we need, to evaluate M a b, is the dependence of the dihedral angles 9 e on the edge lengths. To obtain this 
dependence, we start with the following well known expression of the dihedral angles, 



sm( 



3 e e v 

2A 1 A 2 [ 



(A2) 



where A\ and A 2 are the areas of the two triangles sharing the edge e. All the quantities appearing above can be 
expressed in terms of the (squared) edge lengths, using the Cayley matrix, 



in) 



( 1 1 

' 1 i\ 

1 i\ 

V i e 2 f 2 

\ 1 l n l 2n-l 



C 2 

"o J 



(A3) 
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Indeed, we have 



V 2 

V (n) 



(_!)«+! 



(A4) 



Using this formula for n = 3 to get the volume, and for n — 2 to get the areas entering (| A2|> . we can compute all the 
derivatives of the Regge action. For the equilateral tetrahedron at hand, with £q = Jo + 5 and cos?9 = — i, we have 
the following values: 



/ 



M„ 



6$l 


_V2 
U a 
13 

209 
54^2^ 
1715 



t) 

_V2 
to 
1 

~V2£f, 
33 



_V2 
3 

21 



13 

5 

21 



209 
' 54V5^ 
33 
4\/2^ 



1715 
108\/2# 



V 



(A5) 



These are all the necessary quantities to compute the NLO, as in section ITITl and the NNLO, as in section llVl 

In Section IVIII we generalized the geometric setting and considered an isosceles tetrahedron. In this case, we have 
two different dihedral angles, corresponding to the space edges £q, and the bulk edges It- 



and 



cosflo = 



cos v t 



AC 2 — %f 2 

sin* = 2V2£ n 



Using these, and proceeding as above, the matrix of derivatives of the Rcggc action is 

/ 4e t £ t + 2e £ — s ^- 



(A6) 



(A7) 



dj?dj% 



]o,0t 



170 
V2 



m-z, 



\ 



(A8) 



These are the only relevant terms for computing the LO (|59|1 . 



APPENDIX B: BESSEL FUNCTIONS 



The modified Bessel functions of the first kind are defined by the following integral, 

Ij(z) = - [ d(/)e zcos<l> cos(j(j>). (Bl) 
n Jo 

Using this definition, it is straightforward to verify 1411) and l|42|) . For the calculations of this paper, we have z = 1/a 
are interested in the a 1— > limit, i.e. z 1— ► 00. In this limit, we can evaluate i|Bl|) by a saddle point approximation. 
Firstly, we write the cosine in exponential form, so that (|B1|) becomes a sum of two exponentials, with arguments 
z cos 4>±ij(f). Consider the positive exponential first, z cos <p+ij(j). The stationary point is given by solving the equation 
sin</>o = i-. To do so, we define (f>o = iipo, and use sm(iipo) = 2sinh(?/;o), so that ipo satisfies sinh?/>o = f — "00 m the 
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z i— > oo limit. Consequently, cos^o = y l + (f) 2 — 1 + f(§) 2 - Expanding the exponent of (|Blj) around </>o we thus 
obtain 



j i 

COS</) +«-0O - 77COS0 O (<^- <fo) 2 

z 2 



i-k J -r--( 

2 z 2 



Jn2 



(B2) 



It is easy to check that the negative exponential gives the same result; using it in IBlfl we have 



2- 



■f(0-*i) 2 



e e 



(B3) 



Consequently, in this limit Ij(z) Ij+\(z). 

From these results, it is straightforward to derive (|45() . 
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